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Abstract 

The ground state configurations of the solution to Skyrme's topological soliton model for systems 
with baryon number larger than 1 are well approximated with rational map ansatze, without 
individual baryon coordinates. Here the canonical quantization of the baryon number 2 system, 
which represents the deuteron, is carried out in the rational map approximation. The solution, 
which is described by the six parameters of the chiral group SU(2)xSU(2), is stabilized by the 
quantum corrections. The matter density of the variational quantized solution has the required 
exponential large distance falloff and the quantum numbers of the deuteron. Similarly to the axially 
symmetric semiclassical solution, the radius and quadrupole moment are, however, only about half 
as large as the corresponding empirical values. The quantized deuteron solution is constructed for 
representations of arbitrary dimension of the chiral group. 

PACS numbers: 12.39Dc, 14.20Pt 
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I. INTRODUCTION 



The classical ground state solutions to Skyrme's topological soliton model for baryons 
[lj], with baryon number (B) larger than 1, have intriguing geometrical structure, with 
polyhedral symmetry 2J. The simplest example is the system with B = 2, which has 
axial symmetry in agreement with the description of the deuteron based on a quantum 
mechanical Hamiltonian for the interacting two-nucleon system |4j. Simple rational map 
ansatze, which provide remarkably accurate approximations to the classical ground state 
configurations, have been derived for several systems with baryon number larger than 1 . 
These rational maps represent formal generalizations of Skyrme's hedgehog ansatz for the 
single baryon. Here the rational map ansatz for the B = 2 skyrmion is employed to carry 
out the canonical quantization of the B = 2 solution, which represents the deuteron. 

The Lagrangian density of the Skyrme model is chirally symmetric under constant 
SU(2)xSU(2) transformations. The parameters of this symmetry group are treated as the 
collective dynamical coordinates in the quantization procedure. 

The semiclassical quantization procedure developed in ref. j(| for the nucleon employs 
only that half of the parameters of full chiral symmetry group, which correspond to the 
diagonal subgroup, and therefore it describes only baryons with equal spin and isospin. The 
description by the SU(2) Skyrme model of such states, which have unequal spin and isospin 
as the deuteron, requires six dynamical variables in order to allow consideration of separate 
rotation of the generators of the SU(2) group and the spatial vectors Q|. 

Because the generators are irreducible tensors, these rotations are not independent, how- 
ever, and therefore cannot be not be used for canonical quantization. In the semiclassical 
quantization procedure the classical skyrmion is treated as a rigid body with the conse- 
quence that the quantum mechanical rotation contributes only positive terms to the energy 
functional. There is then no stable variational solution to the quantized Hamiltonian. To 
obtain a stable variational quantized solution one may draw on the canonical quantization 
procedure, which has been developed by K. Fujii et al. [a] for the SU(2) Skyrme model. The 
treatment of the dynamical field variables of the Skyrme Lagrangian density as quantum- 
mechanical variables ab initio generates negative quantum corrections and also stable quan- 
tum solitons |9| . These energy of these quantum solitons depends on the dimension of the 
representation of the symmetry group in contrast to the semiclassical case lfll 
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The canonical quantization procedure developed below employs the two sets of three 
Euler angles, which correspond to left and right chiral rotation groups as the six collective 
coordinates. The resulting canonical angular momentum operators lead to compact forms 
for both the Lagrangian and Hamiltonian of the biskyrmion. The two sets of independent 
angular momentum operators allow construction of the eigenstate of the biskyrmion from 
the eigenstates of two subsystems. For the deuteron the subsystems are the neutron and 
proton, which form the state with common spin S = 1 and isospin T = 0. The approach 
generalizes to dibaryons, which may be constructed from neutrons and protons as well as 
from A resonances. 

The matter density of the quantum soliton falls off exponentially at long range in contrast 
to power law falloff of the classical solution. The inverse of the length scale of this exponential 



falloff for the B = 1 skyrmion corresponds to the pion mass 



D 



In the case of the deuteron it 



should correspond to 2^m^B, where B is the binding energy and the nucleon mass. It 
is shown here numerically in the rational map approximation that for the variational ground 
state the matter density falls off at roughly this rate, as required. 

The approximate quantum soliton for the deuteron derived here describes the rotational 
quantum corrections appropriately, but not the large distance solution of two well separated 
single skyrmions. This is revealed by the magnitude of its radius and quadrupole moment, 
which are only about half as large as the corresponding empirical values. The semiclassical 
solution shares these features 0. 

The present manuscript is organized as follows. In Section 2 the classical rational map 
ansatz for soliton with baryon number 2 is generalized to representations of arbitrary di- 
mension. In Section 3 this soliton (biskyrmion) is canonically quantized with six collective 
variables, which correspond to the parameters of chiral symmetry SU(2)xSU(2) group. The 
expressions for the electric form factor, quadrupole moment and rms radius of the deuteron 
are presented in Section 4. The numerical results for deuteron observables are discussed in 
Section 5. 

II. THE CLASSICAL AXISYMMETRIC SOLITON 

The Skyrme model is a Lagrangian density for a unitary field U(x, t), which may described 
by any representation of the SU(2) group. In a general reducible representation the most 
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compact expression of the unitary field £/(x, t) is as a direct sum of Wigner's D matrices for 
the irreducible representations of arbitrary integer or half integer dimension j [lol |: 

U(x,t)=Y / (BD i (a(x,t)). (1) 

3 

The D matrices depend on three unconstrained Euler angles ol = (a 1 , a 2 , a 3 ). 
The chirally symmetric Lagrangian density has the form 

C[U(x,t)\ = -&Tt{R^} + y^Tr{[^,^] 2 }, (2) 

where the "right" current is defined as 

= {d^U)U\ (3) 

and f n (the pion decay constant) and e are parameters. 

The static variational solutions to classical Skyrme model for baryon number B = 2 
(biskyrmion) have been derived numerically in refs. I, Q. In ref.fl the following simple 
rational map ansatz, which preserves the axial symmetry of ground state solution for the 
biskyrmion, was found to give an approximation to the ground state energy, with an accuracy 
of better than 3 per cent: 



c 



Here J a are SU(2) generators, which may be defined in representations of arbitrary dimen- 
sion. The scalar function F R {r) is the "chiral angle" for the biskyrmion, which is determined 
by the variational equation of motion. The circular components of the unit vector n are 

~-i sin 2 $ 2i 

n+i = -n = - e il¥ \ 

^(l + cos 2 ^) 

.0 2 cos?? , , 

1 + COS^ V 

-+i sin 2 $ _ 9iu , 

n_i = -n +l = —= e v . 

^(l + cos 2 ^) 

Substitution of the ansatz (£Q) in the Lagrangian density (J2J) yields the following expression 
for the mass density of the classical biskyrmion: 



M d = ^{f 2 n [F' 2 (r) + 21 sin 2 F R (r 



N 
I 

21 1 i 

+ — sin 2 F R {r) [F l2 {r) + - sin 2 F R (r)]j. (6) 



Here X is the Gaussian curvature, defined as 

1= 4sinH , (7) 
r 2 (l + cos 2 tf) 2 1 ; 

In contrast to the hedgehog ansatz for B = 1 this mass density depends on both the polar 

angle d and the radius r. The classical Lagrangian density depends on representation only 

through the overall factor N = | YljjU + 1) (2jf + 1), where 2j + 1 is the dimension of the 

representation, and which may be absorbed by renormalization of the model parameters 

0. 

The requirement that the soliton mass be stationary yields the following differential 
equation for the chiral angle: 

f r(t) (f 2 + 4 sin 2 F R (r)) + 2F' R 2 (f) sin 2F R (f) + 2rF' R (r) 

- sin 2F R (f) (2r 2 + (- + tt) sin 2 F R (r)) = 0. (8) 

Here the dimensionless variable f is defined as f = ef n r. At large distances r —>■ oo this 
equation reduces to the asymptotic form 

r 2 F R (r) + 2fF' R (f) - AF R (f) = 0. (9) 

The solution to the asymptotic equation Q falls of with an algebraic power of distance: 

F R (r) ^Cf- 1 ^. (10) 

The falloff rate is somewhat larger for the biskyrmion than for the hedgehog ansatz for 
B = 1, as the power of f in Eq. ()1(JJ) is —2.56 , whereas in the case of B = 1 it is —2. After 
renormalization by the factor N, the biskyrmion baryon density takes the form 

B°(r) = -^F R (r) sin 2 F R (r). (11) 

III. CANONICAL QUANTIZATION WITH SIX COLLECTIVE VARIABLES 

The Skyrme Lagrangian (j2J) is symmetric under chiral SU(2)xSU(2) transformations. 
The canonical quantization of the classical soliton solution (jHj) can be achieved by means of 
collective coordinates that separate the time dependent variables from those that depend on 
the spatial coordinates: 

U(x, a(t),(3(t)) = A (a(t)) U R (n, F R {r))B^ ((3(t)) . (12) 



Here the two sets of three Euler angles a(t) = a 2 (t), a 3 (t)) and (3(t) = 

(t) , /? 2 (t) , /3 3 (t)) are those for the two SU(2) groups respectively. In the canonical quan- 
tization the Skyrme model is considered quantum mechanically ab initio. The collective 
coordinates ot(t),/3(t) and velocities ct(t),/3(t) are treated as dynamical variables with the 
commutation relations 



[a k ,a l ] = -i R f kl (cx,cx), [a\(3 l ] = -i R f kl (a, (3), 



(13) 



i*A/3,/3), 



(3 k , a 1 



The functions R f kl are defined in (|2*H> below. 

Substitution of ansatz ([T2|) into the Lagrangian density (J2J yields the quantum La- 
grangian, which is quadratic in the generalized velocities: 



L =^a k g k i{ot, a)a l + ^a k g kl {a, /3)$ l + ^$ k g M {f3, a)a l 
+ ^ k g kl ((3,f3)$ l + [(«,/3)°-order term]. 



(14) 
(15) 



Here the coefficients g k i are defined as 



g H (a,a) = ci a \oc) l R ab {F R )C?\ct), 
g u (P,fl = C' k {a \f3) 1 R ab (F R )C' l {b \f3), 
g]a (a,p)=g lk (J3,a) = Ci a \ a ) 2 R ab (F R )Cl ib \l3). 



(16) 



The C' k s and their inverses C'A are functions of the dynamical variables, which appear in 
the differentiation of the Wigner D matrices [lo| : 



d 



1 



<t\") Dl m ,{cc)(jrn'\j a \jn). 



The matrices i^Rab{F R ) are antidiagonal 



(17) 



iR ab {F R ) = (-l) a 1R{F R )5 a ^ 



lRab\F R ) = (— 1)° %R(F R )5 a -b, 



(18) 
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and have the matrix elements 



lR(F R ) = -Y-^j~J drT 2 (1(4 - vr)^ 2 + ± sin 2 F R ) , (19) 
fR(F R ) = Y -J^J- J d ^ 2 sin 2 F R + l n cos2F R Fg + sin 2 2F R 
° 2 R(F R ) = Y — ^j— y drT 2 ((4 - vr) sin 2 Fr — (4 — vr) cos 2F R F' R 2 

These matrix elements contain the infinite integral 



Y= ^ki df?2 ' (20) 

which drops out from the generalized moments of inertia and the mass density, when taken 
to infinite, but which is convenient to retain formally in the intermediate steps. 

The infinite terms arises in the quadratic term in the Lagrangian density, when the left 
and right rotations are unequal: 

Tr R°R « TiAA ] AA ] + Tr BB^BB ] 

-Tr AA^U Q BB ] Ul - Tr UqAA^U^BBI (21) 

Here only the terms, which contain 6t or (3, and which are important for commutation 
relations are considered. The infinite terms in (|T9*|) arise from the terms on the r.h.s. of (|2ip. 
which are independent of the spatial coordinates. In the case when A — B, Uq — *■ 1 when 
r — > oo the infinities disappear from the Lagrangian ()14j) . 

The Lagrangian (fT^j) may be used to define the following canonical momentum operators, 
which are conjugate to the collective coordinates: 

7T fe («) = ^ = ^{a\gu(a,a)} + ^{/3 l ,g kl (a,/3!)}, (22) 
7r*(/3) = ^ = -{/5 / , fe (A/3)} + 2(^'^(A«)}. 

Here the curly brackets denote anticommutators. The canonical commutation relations 

[7T k (cc),a l ] = -i5 k , h [n k ((3),/3 l ] = -i5 k , h [n k (a),p l ] = [n k ((3),a l ] = (23) 
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lead to the system of linear equations for the functions R f kl in (|13p. the solution of which 
can be written in the form 

R f hl (a,a) = C5 ) (a) 1 F ab (F i? )Cj ) (a), 

R f kl (P,f3) = Cfa^F^F^CfoiP), (24) 
*/ w (a,/9) = R f\(3^) = C' k a) {(x) 2 F a \F B )C'l b) {(3). 

Here the antidiagonal matrices 



F ab (F R ) = (-IT ?F(F fl )5 a ,_ 6 , (25) 



2 F ab (F R ) = (-l) a a 2 F(F R )8 a ^ b , 
have the following matrix elements, which in the limit Y — > oo become finite: 

lim fF = lim 

Y—*oo Y— >oo 

lim g F — lim 



x— +oo 



lim ®F = lim 



fR 


2 






-fR 


2 


(fRy-(±Ry- 


di 


° 2 R 


2 


my - my 


a ' 


-1R 


2 


my - my 


Oo' 



\imjF = lim 7W1 ^- 7W1 ^ = --• (20) 
The quantities 

a °=^=ikr dr> ~ 2 sin2 fr{{a fs) + 8 ^ } ' (27a) 

~ /*00 • 2 771 

^ = wrwJo d - 2sin2 ^H i +^)+ 16 ^)' ( 27b ) 

define two different soliton moments of inertia, as appropriate for an axially symmetric 
system. It is convenient to introduce the following angular momentum operators on on the 
hypersphere S 3 , which is the group manifold of SU(2): 

J» = --^{^(aj.q^o)}, (28) 

the components of which satisfy the standard commutation relations and 

[j>),j>(/3)]=0. (29) 
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The coefficients of the quantized Lagrangian (|14j) contains the infinite integrals Y. After re- 
placement of the velocities by the natural angular momentum operators ()29|) . the Lagrangian 
density can be reexpressed as a sum of the angular momentum operators J' {ex.) + J'{f3). By 
means of some lengthy manipulation the Lagrangian density takes the following form in 
terms of these: 

£(«, A h, F R {r)) = -M d + sin 2 F R {l + + Jsin 2 F R ) 

x (l(j'(a) + J'{/3)) 2 + (1 - 1) (4(a) + 4(/3)) 2 (30) 

- "J-(J'(a) + J'{(3)) ■ n) + (1 - (J^(a) + 4(/3))n l ') + AM. 

The last term on the r.h.s., AAi, is the quantum correction to classical mass density, which 
appears on account of the commutation relation (|23J1 . This has the expression 



where 



AM = f n e |_ + _ + _ + d(_ + _ + _U, (31) 



Q 3 = lsm 2 F R {l-hl){l + Fg+lsm 2 F R ), 

o 

Q 4 = i sin 2 F R {l-hl){l + F%), 

Q 5 = isin 2 F iJ [(l + 3n 2 )(l + F^ 2 ) + (l+r i 2 )Jsin 2 F i? ], 

Q 6 = ^ sin 2 F R {1 - h 2 ) [{1 - n 2 )(3 sin 2 F R + 3Fg - 2F' 2 sin 2 F R ) (32) 

+ (l + 3n 2 )Xsin 2 F iJ ], 
Q 7 = ^sm 2 F R {l-h 2 )[{l + 3nl){sm 2 F R + Fg) 

- sin 2 F R {2{1 + hl)F'i + 3Jn 2 )] , 
Qs = ^sin 2 F fi [(3 + 2n 2 + 3^)(sin 2 F fi + F^ 2 )-sin 2 F fi (2(l + n 2 )F j [ 2 

-(1 + An 2 - 3n$)X)] . 

Here only the quantum correction depends on the dimension of the representation of the 
chiral field U through the explicit factor 

, 3 E J j(j + l)(2j + l)(2j-l)(2j + 3) 

2 E,i(i + l)(2i + l) • (33) 

Traditionally the Skyrme model is formulated in the fundamental representation, in which 
j — | and d — 0. 



The operators in the Hamiltonian of the biskyrmion system also are the sum of two 
independent angular momentum operators J' (ex.) + J'(/3)- The terms with the operators 
(j'(a)) 2 , (J'((3)) 2 or j'(ex)J'({3) drop out from the Hamiltonian as they contain coefficients 
with the infinite factor Y in their denominators. The angular momentum operators are 
natural operators for Skyrme model and in terms of them the Hamiltonian operator for the 
biskyrmion becomes: 



where 



H(cx,f3,F R )= Md + ^-(f(cx) + f(f3)) 2 

+ (4-4)(^(«) + J^)) 2 + AM, 



AM 



(34) 



d 3 rAM, (35) 

is the quantum correction to soliton mass. The Hamiltonian is similar to semiclassically 
quantized Hamiltonian of a rotator, with exception for the quantum correction. The normal- 
ized eigenstate vectors for the Hamilton operator (}3~l*j) can be constructed from eigenstates 
of two subsystems with common spin S and isospin T as: 
S T 



h k 

m s m t 



v /(2/ 1 + l)(2/ 2 + l) 



16tt 2 



x 



mi ,7712 ,77k , mo 



h k s 

m\ rri2 m s 



h k T 
m[ m' 2 m t 



D 1 , (a)D 

mimj V / m2m 2 



,((3)\0). 



(36) 



The operators J f (oc) and J' ((3) are "right rotation" operators for the Wigner matrices 
D ll (cx) and D £2 (f3). The biskyrmion with different S and T can now be constructed from 
states with the quantum numbers of the nucleons and the A resonances. The eigenvalue of 
Hamiltonian operator gives the mass of quantum biskyrmion as 



M d = M d + AM + - 



— T(T+ 1) + (—- — )m 
a% Vao ai 



(37) 



which depends only on isospin T and isospin projection m t . For the deuteron T = 0, and the 
variation of the mass ()37|) gives the integro-differential equation for the chiral angle Fr(t) 

SM d 
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0. 



(38) 



This explicit expression for the deuteron is identical to the corresponding equation for 
dibaryons At large distances the equation (j^S} reduces to the asymptotic form 



^2 jpn 



+ 2rF' R - (4 + m 2 ? 2 ) F R = 0, 



(39) 



where 



and 



m 2 



A -it 



K + 7)-^I (r(T + l)-m? + 2)7r + 2 



/V 3 



4 y 45f 

4 — 7T 7T 



7. 7.2 



+d 



2(4 -?r) ~ 



^3 



Q 6 + 



4 - 7T 



+ 



7T 



7T ~ 
7T ~ 



Q7 + "^jQs 



d 3 rQ k . 



4 - 7T 

4aoSi 



(40) 



(41) 



The factor m describes the falloff rate of the chiral angle at large distances: 



Fr (r) 



C[— + — e 



(42) 



The related quantity m = ef n m describes the asymptotic falloff exp(— 2mr) of biskyrmion 
mass density like Yukawa pion cloud for nucleon. 



IV. STRUCTURE OF THE QUANTIZED DEUTERON SOLUTION 

The electric and quadrupole form factors Gc(Q 2 ) and Gq(Q 2 ) of the deuteron state 
are obtained as the matrix elements of the spin-scalar and spin-tensor parts of the time 
component of the electromagnetic current operator. For the isospin deuteron state this 
is given by the anomalous baryon current. The matrix element is evaluated between the 
deuteron states in the Breit frame, which is defined by p + p' = (Q): 

(d ,m' s p'\j°{r = 0)\dm s p)= Gc(Q 2 )5 mam > 

+ 6M! G QW 2 )^>(3?V - Q 2 8 ab )Ut ms . (43) 

Here q = p' — p is the momentum transfer, Q 2 = —q 2 , is the mass of the deuteron, and 
U ms a is the unitary matrix that relates the Cartesian and spherical bases. 
The expression for the electric form factor is: 

G C {Q 2 ) = \J d 3 r Jo (gr)i3°(r), (44) 
11 



where jfc(gr) is the spherical Bessel function of k-th order. The quadrupole form factor is 
correspondingly 

s /if 2 r 

°Q(Q 2 ) = 2~fj d3r(1 - 3cos2 ^(^)i3°(r). (45) 
The mean square charge radius and the quadrupole moment are defined as 

r 2 = -6^G C (Q 2 ), Q = M 2 d G Q (0). (46) 

It follows that 

(r 2 ) ch = \J d 3 rr 2 B°(r) (47) 

and that 

Qd = Yoj dhr2{1 " 3cos ^) i3 °( r )- ( 48 ) 

The matter radius of the deuteron solution, (r 2 ) m , may in turn be determined from deuteron 
mass distribution as 

(r*)m = Tr I d 3 vr 2 (M(r) + AM(r)). (49) 
Md J 

V. NUMERICAL RESULTS AND DISCUSSION 

The numerical value for the rate m (|4())l . at which the mass density of the solution 
decays with distance provides the key test of the phenomenological gain in the canonically 
quantization of the ground state solution with the quantum numbers of the deuteron. In 
order to match the falloff rate of the matter density that corresponds to the deuteron wave 
function, it should equal 2yfBM n , where B is the binding energy and M n the nucleon mass. 
The empirical value of this quantity is 91.4 MeV. This test requires numerical solution of the 
variational problem for the quantized deuteron state. For this purpose the two parameters 
f n and e of the Lagrangian density of the Skyrme model have, however, to be determined 
first by fits to two empirical nucleon observables. 

The procedure adopted here was to first determine these two parameters by using the 

chiral angle of the classical Skyrme model, which is independent of both model parameters 

rn 

and the representation [10], so that the empirical values of the mass (M n = 939 MeV) and 
isoscalar radius (0.77 fm) of the nucleon are reproduced. These parameters were then used 
in a numerical solution of integrodifferential equation for the chiral angle [9j, which does 
depend on the dimension of the representation in the case of the quantized skyrmion. That 
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solution was subsequently used to determine new values of f w and e. This procedure was 
iterated until a convergent solution was obtained. 

The values for the parameters found by this method are listed in Table for a set of 
representations of the SU(2) group with different dimension j. The chiral angle for the 
deuteron state, Fr(t) was then determined by self consistent numerical variation of the 
energy expression ([37)1 . for the representations listed in Table H] 

The numerical results for the deuteron mass Md, the binding energy AE = Md — 2M n , 
matter radius r m , charge radius r ch , electric quadrupole moment Qd and the mass rh pOjl. 
which describes the falloff rate of the deuteron mass distribution at large distances, are listed 
in Table for the irreducible representations with dimension j = 1/2, 1 and 3/2 as well as 
for the reducible representation 1 © 1/2 © 1/2. The quantum correction to the ground state 
energy is similar in size to that found in ref. 3| by an entirely different approach. The value 
of the falloff mass rh, is in all cases considered of the same order of magnitude as the quantum 



mechanical value 2y/BM n , and in the case of the three dimensional representation actually 
agrees with that value. The value is also close in the case of the reducible representation. 
It is interesting that these two representations are also those, which give the best values for 
the falloff rate for the matter density of the nucleon, which should be of the order of the 
pion mass 

The shape of deuteron is represented by the mass density distribution M. d + AM.. The 
equidensity surface displayed in Fig. ^ is roughly toroidal. The maximum value of mass 
density is 1150 MeV-fm -3 as shown in Fig. 121 The density maxima form a ring with a 
diameter of 1.424 fm. The baryon density distribution has a similar shape. An analogous 
toroidal structure in the deuteron has been shown to arise in the quantum mechanical 
treatment of the two- nucleon system with a realistic interaction Hamiltonian [l^ . 

The calculated nonrelativistic charge form factor is shown in Fig.|3]for the representations 
with dimension j = |, 1 and |. The calculated form factor has the same qualitative features 
as the empirical form factor value, although the charge radii are much too small. The 
calculated charge and matter radii as well as the quadrupole moments are listed in Table HI 
These values are only about half as large as the corresponding empirical values, a result 
which is similar to that found for the semiclassical axisymmetric skyrmion description of the 
deuteron. More realistic values for these static parameters, which represent the large scale 
features of the deuteron obtain with Skyrme's product ansatz for the two-nucleon system 
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TABLE I: The predicted static deuteron observables in different representations with fixed empirical 
values for nucleon isoscalar radius 0.77 fm. and mass 939 MeV. 



3 


1/2 


1 


3/2 




Expt. 


u 


57.68 


56.53 


55.71 


56.82 


93 MeV 


e 


4.325 


4.08 


3.79 


4.13 




M d 


1868 


1926 


1998 


1907 


1876 MeV 


AE 


-10 


48 


120 


29 


-2.22 MeV 


r c h 


1.10 


1.17 


1.24 


1.15 


2.13 fm 




1.22 


1.29 


1.35 


1.27 


1.97 fm 


Q d 


0.140 


0.158 


0.177 


0.152 


0.286 fm 2 


m 


54.6 


90.5 


110.0 


82.3 


91.4 MeV 



0. 
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FIG. 1: Equidensity surface for the quantized deuteron 




FIG. 2: Mass density distribution of the quantized deuteron 



16 



10* • 



-1 — I — 1 — I — . — I — . — I — . — I — . — I — . — I — . — I — 1 — I 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 

Q 2 (GeV/c) 2 



IG. 3: Electric form factor of the quantized deuteron solution. The experimental data are from 
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